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We compute the energy density radiated by a quark undergoing circular motion in strongly coupled 
M = 4 supersymmetric Yang-Mills plasma. If it were in vacuum, this quark would radiate a beam of 
strongly coupled radiation whose angular distribution has been characterized and is very similar to 
that of synchrotron radiation produced by an electron in circular motion in electrodynamics. Here, 
we watch this beam of gluons getting quenched by the strongly coupled plasma. We find that a 
beam of gluons of momenta ~ q 3> ttT is attenuated rapidly, over a distance ~ q 1 ^ 3 (nT)~ 4 ^ 3 in a 
plasma with temperature T. As the beam propagates through the plasma at the speed of light, it 
sheds trailing sound waves with momenta < nT. Presumably these sound waves would thermalize 
in the plasma if they were not hit soon after their production by the next pulse of gluons from the 
lighthouse-like rotating quark. At larger and larger q, the trailing sound wave becomes less and less 
prominent. The outward going beam of gluon radiation itself shows no tendency to spread in angle or 
to shift toward larger wavelengths, even as it is completely attenuated. In this regard, the behavior 
of the beam of gluons that we analyze is reminiscent of the behavior of jets produced in heavy ion 
collisions at the LHC that lose a significant fraction of their energy without appreciable change in 
their angular distribution or their momentum distribution as they plow through the strongly coupled 
quark-gluon plasma produced in these collisions. 



I. CONTEXT 



Jet quenching is one of the most striking phenomena 
seen in heavy ion collisions at the Relativistic Heavy 
Ion Collider (RHIC) and at the Large Hadron Collider 
(LHC). The basic picture behind a suite of observables 
measured in these collisions is that energetic partons pro- 
duced via rare hard-scattering processes in heavy ion col- 
lisions lose a significant fraction of their energy as they 
and the spray of partons they fragment into plow through 
the strongly coupled plasma produced in the same col- 
lisions. At the LHC, jet energies are high enough that 
the jets can be detected calorimetrically event-by-event, 
and the phenomenon of jet quenching is manifest in sin- 
gle events with, say, a jet with an energy greater than 
200 GeV back-to-back with a jet with an energy less 
than 100 GeV [THS]. (It is improbable that a pair of 
jets will be produced such that each travels the same dis- 
tance through the plasma and loses the same amount of 
energy, so back-to-back pairs of jets with unbalanced en- 
ergies are the norm.) At the LHC, the attenuated jets 
(i.e. the lower energy jets in the unbalanced pairs) look 
remarkably like jets produced in vacuum, with angular 
distributions and momentum distributions (i.e. fragmen- 
tation functions) that are to date indistinguishable from 
those of jets with the same energy produced in proton- 
proton collisions [H HI [5] or in peripheral heavy ion colli- 
sions [5]. The energy lost from the jets emerges instead 
as an excess of soft particles (momenta < 1 GeV [2]) at 
large angles (> 45° [5]) relative to the jet direction. 
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The observation of jets that have lost a large frac- 
tion of their energy was no surprise; it was surprising, 
however, to see the attenuated jets emerging without 
any visible softening of their fragmentation functions and 
without any visible broadening of their angular distribu- 
tion. In the limit of high parton energy E, the domi- 
nant energy loss process for an energetic parton plowing 
through quark-gluon plasma with temperature T is gluon 
bremsstrahlung [7]-[9] , radiating gluons with energy ui and 
momentum transverse to the jet direction k± that satisfy 
E >• lj 3> kj_ 3> ttT [BHT2"]. This set of approximations 
underlies all analytic calculations of radiative energy loss 
to date. The (perhaps naive) expectation based upon 
these considerations is that at least some of the energy 
lost by the high energy parton should emerge as relatively 
hard particles (since u) 3>3> 7rT) near the jet direction 
(since u) 3> k±), resulting in a jet whose angular dis- 
tribution has been broadened and whose fragmentation 
function has been softened. The data from the LHC have 
been stimulating more sophisticated implementations of 
these considerations [I3ll21j . but it is also possible that 
the partons produced in LHC collisions are simply not en- 
ergetic enough for this picture to apply at all — since it 
is based upon the premise that the QCD coupling evalu- 
ated at the scale k± (which, recall, is <C^C E but ^> ttT) 
is weak even if the physics at scales ~ ttT is strongly 
coupled. Given that this separation of scales may not 
be applicable, it behooves us to analyze jet quenching 
or models of jet quenching in strongly coupled plasmas 
in contexts where reliable analyses are possible. Even 
if these analyses yield only qualitative insight, they can 
be useful as benchmarks and as guides to how to think 
about the physics. 

The simplest strongly coupled plasma that we know 
of is that found at nonzero temperature in strongly cou- 
pled Af = 4 supersymmetric Yang-Mills (SYM) theory 
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in the large number of colors (N c — »• oo) limit, which 
is dual to the five-dimensional anti-deSitter black hole 
(AdS-BH) [2"2"H2"4"] . (For a recent review of the by now 
many ways in which this theory and its many cousins 
have been used to gain insight into properties of strongly 
coupled plasma and phenomena in heavy ion collisions, 
see Ref. [IS].) Unfortunately, it is not possible to liter- 
ally study jet quenching in Af = 4 SYM theory because 
hard scattering in this theory does not produce jets [2"o] . 
Nevertheless, many authors have used the strongly cou- 
pled plasma of Af — 4 SYM theory to gain relevant in- 
sights, for example by studying the energy loss and mo- 
mentum diffusion of a heavy quark plowing through the 
plasma [27H3T] as well as the wake it produces [32H39] . 
the stopping distance of an energetic light quark or gluon 
in this plasma |4T)H4"2"] , and the value of the jet quenching 
parameter [43] [44] , the property of the strongly coupled 
plasma that enters into the calculation of jet quenching if 
the scale-separation of the previous paragraph does turn 
out to be valid. However, none of these calculations yield 
insights into how an attenuated jet can emerge with- 
out any visible softening of its fragmentation function 
or spreading of its angular distribution. The one that 
comes closest is the calculation of the stopping distance 
of a light quark [40H42]. 



Note that it does not make sense to interpret the data 
as saying that a single quark loses a lot of energy in the 
medium and then emerges in isolation and fragments into 
an ordinary-looking jet, since what would emerge is a 
nearly onshell quark, which would not fragment into a jet. 
If we are to gain insight into jets in heavy ion collisions 
from a strongly coupled perspective, we should imagine 
that the initial hard parton fragments quickly into a pro- 
tojet of some sort and then ask how this protojet in- 
teracts with, and loses energy in, the strongly coupled 
plasma. Although there is a sense in which the analyses 
of Refs. [4"0H4"2"] provide answers to this question, these re- 
sults are sensitive to details of the initial conditions, and 
we anyway do not expect strong coupling methods to de- 
scribe the initial fragmentation of a hard parton into a 
protojet. Instead, we shall assume that this initial stage 
proceeds conventionally, as described successfully in per- 
turbative QCD, and then ask how the protojet interacts 
with, and loses energy in, the strongly coupled plasma. 
Recent work |45j offers a new way to gain a perspec- 
tive on this question, since it provides a way to produce 
a beam of gluons whose angular distribution and wave- 
lengths are well understood in vacuum. In this paper, we 
shall shine such a beam of gluons through the strongly 
coupled plasma at nonzero temperature and watch this 
beam rapidly get attenuated — without any apparent 
broadening of its angular distribution or lengthening of 
its wavelength! 



II. INTRODUCTION 

The trick by which a beam of gluons can be produced 
in Af = 4 SYM theory is to consider a test quark un- 
dergoing circular motion with radius Rq and velocity ft 
(and hence angular velocity Q = f>Ro) in the vacuum 
of this theory [J5J HB]. At both weak coupling (where 
the calculation is done conventionally) and strong cou- 
pling (where the calculation is done via gauge/gravity 
duality) the radiation that results is remarkably similar 
to the synchrotron radiation of classical electrodynamics, 
produced by an electron in circular motion |45) . In par- 
ticular, as the limit of ultrarelativistic motion is taken 
(7 —> 00 where 7=1/ yl — j3 2 ) the lighthouse-like beam 
of radiation becomes more and more tightly collimatcd 
in angle (it is focused in a cone of angular extent ~ 1/7) 
and is composed of gluons and scalars with shorter and 
shorter wavelengths (the pulse of gluons in the beam has 
a width ~ Rq/j 3 in the radial direction in which it is 
moving). The emitted radiation was found to propagate 
outward at the speed of light forever without broadening 
either in angle or in pulse width, just as in classical elec- 
trodynamics [4"5H5"T] . At weak coupling, the slight differ- 
ences in the angular distribution of the power radiated to 
infinity relative to that in classical electrodynamics can 
be attributed to the fact that scalars are radiated as well 
as gluons [35] ■ And, at strong coupling the angular dis- 
tribution is identical to that at weak coupling 45, 48, 50 . 

We shall compute the energy density radiated by a 
quark undergoing circular motion in strongly coupled 
Af = 4 SYM theory at nonzero temperature, allowing 
us to watch what happens as the lighthouse beam of glu- 
ons and scalars is attenuated as it shines through the 
strongly coupled plasma. We shall gain analytic under- 
standing of the length scale over which the energy of the 
beam is attenuated by the plasma. And, by inspection of 
the energy density, which we obtain numerically, we shall 
see that as the beam is attenuated it does not broaden 
in angle or redden in wavelength. 

The rate at which a quark undergoing circular motion 
through the plasma of strongly coupled Af = 4 SYM 
theory loses energy was studied previously in Ref. |52j . In 
this analysis, two distinct regimes were found, depending 
on whether 
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is 3> 1 or -C 1. For E>1, the energy loss rate is given 
by the generalized Larmor formula 



dE 
~dt 



rati 



2ir ' 



(2) 



where is the quark's proper acceleration and A = g 2 N c 
(with g the gauge coupling) is the 't Hooft coupling which 
can be chosen at will since in this conformal theory it 
does not run and which we take as large since we wish to 



3 



study strongly coupled plasma. It was shown many years 
ago by Mikhailov that the energy loss rate of a quark in 
circular motion in the vacuum of strongly coupled Af = 4 
SYM theory is given by ([2| [S3] and so in the strongly 
coupled plasma at T ^ 0, in the H> 1 regime we expect 
to see the radiation of beam of synchrotron-like radiation 
as in vacuum |45) , and the subsequent attenuation of this 
beam. When 3 < 1, on the other hand, the energy 
loss rate is that due to the drag force exerted by the 
strongly coupled hydrodynamic fluid on a quark moving 
in a straight line with velocity j3 [52] , namely [57] [55] 



dE 
~dt 



drag 



VA 
2^ 



(ttT) 2 /3 2 7 



(3) 



Notice that the parameter 3 which governs which expres- 
sion for the energy loss rate is valid is simply the ratio of 
the rates appearing in Eqs. ^ and pb. In this respect 
it is as if both hydrodynamic drag and Larmor radiation 
are in play with the larger of the two effects dominating 
the energy loss, but this simplified picture is not quanti- 
tatively correct because where 5 ~ 1 the energy loss rate 
is less than the sum of Eqs. ^ and ^ [52 . Although 
our principal interest is in the 5 > 1 regime, where we 
can study the quenching of a beam of synchrotron glu- 
ons, it will prove instructive to look at 5 < 1 and S ~ 1 
also as in these regimes the hydrodynamic response of 
the plasma — i.e the production of sound waves — is 
more readily apparent. 

Unlike in vacuum, in the plasma at nonzero tempera- 
ture the energy disturbance created by the rotating quark 
can excite two qualitatively distinct modes in the energy 
density; a sound mode which at long wavelengths travels 
at speed c s = 1/ y/3, and a light- like mode which propa- 
gates at the speed of light. The relative amplitude of each 
mode depends on the trajectory of the quark. We find a 
correlation between the energy loss mechanism and the 
relative amplitude of the light-like and sound modes that 
can be anticipated from the results of Ref. [52]. When 
5 < 1 the dominant modes that are excited are sound 
waves. When S > 1 the dominant modes that are ex- 
cited propagate at the speed of light. Interestingly, by 
studying the S ~ I regime we shall see that as the pulse 
of radiation moving at the speed of light is attenuated in 
energy, it sheds a sound wave. 

We shall use gauge/gravity duality to calculate the en- 
ergy density that results when a test quark is moved 
through the strongly coupled plasma in a circle with 
E < 1 or H ~ 1 or 5 > 1. The gravitational dual of 
the undisturbed plasma is the (4+I)-dimensional AdS- 
Schwarzschild (AdS-BH) geometry. Setting the AdS cur- 
vature radius L = 1, we choose coordinates such that the 
metric of the AdS-BH geometry is 



ds 2 = 



[-f(u) dv 2 + dx 2 - 2 dv du] , 



(4) 



where f(u) = I — (u/uh) 4 . These coordinates are general- 
ized infalling Eddington-Finkelstein coordinates; lines of 



constant time v represent infalling radial null geodesies. 
The event horizon of the geometry is located at u = u^, 
with T = (iruh) -1 the temperature of the equilibrium 
strongly coupled Af = 4 SYM plasma that this metric 
describes. We shall do the calculation in units where 
ttT = I, meaning Uh = 1- (However, in order to fa- 
cilitate comparison to results at T — we shall report 
results instead in units where Rq — I.) The boundary of 
the AdS-BH metric, corresponding via the holographic 
correspondence to the ultraviolet limit in the boundary 
Af = 4 SYM quantum field theory, is at u — > 0. We begin 
in Section III by reviewing the calculation of the shape of 
the string that hangs "down" from the rotating quark at 
the ultraviolet boundary of the AdS-BH toward its hori- 
zon, spiralling around and around infinitely many times 
just above the horizon. This string profile was obtained 
in Ref. 52J , but we rederive it in the coordinates that we 
shall use in Section IV. In Section IV, we solve the bulk- 
to-boundary problem, finding the energy density in the 
boundary Af — A SYM theory plasma that this spiralling 
string describes. We describe our results in Section V 
and in Section VI we present analytic arguments that 
allow us to understand all of their qualitative features, 
and in Section VII we return to the context with which 
we begun, marvel at the qualitative resemblance between 
our results and jet quenching in heavy ion collisions at 
the LHC, and speculate on how this resemblance could 
be made more quantitative. 

There are preliminary indications that the energy lost 
by the lower energy jets (composed of longer wavelength 
gluons) studied at RHIC ends up in soft hadrons mov- 
ing in directions correlated with the initial parton's path 
making the jets appear broadened in angle [54H58] . rather 
than far outside the jet cone as at the LHC [2]. We close 
by observing that our results illustrate a natural way for 
such a distinction between jet quenching at RHIC and 
the LHC to arise, since in our steady-state calculation 
we find that our beam of gluons excites less of a sound 
wave trailing behind the beam pulse — aka soft particles 
going in roughly the jet direction — when the beam is 
composed of shorter wavelength gluons than when it is 
composed of longer wavelength gluons. We discuss this 
in Section VI. 



III. STRING DYNAMICS 

The dynamics of classical strings are governed by the 
Nambu-Goto action Sng — —To J dr do~y/—g , where 
To = \/\/(2ttL 2 ) is the string tension, a and r are the 
worldsheet coordinates of the string, and g = det g a b 
where g a b is the induced worldsheet metric. The string 
profile is determined by a set of embedding functions 
X m {t,<t) that specify where in the spacetime described 
by the metric ^ the point (r, a) on the string world- 
sheet is located. The induced world sheet metric is given 
in terms of these functions by g a b — d a X ■ di,X, where 
a and b each run over (t, a) . For the determinant we 
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obtain 



- g = (d T X ■ d a Xf - {d T X) 2 (d a Xf 



(5) 



We choose worldsheet coordinates r = v and a = u. 
As we are interested in quarks which rotate at constant 
frequency about the z axis, we parametrize the string 
embedding functions via 



X [v, u) = (v,r s (v, u),u) 



(6) 



where in spherical coordinates {r, 9, ip} the three-vector 
r s is given by 



r s (v,u) = (R(u),%,ip(u)+Q.v) 



(7) 



With this parameterization, the Nambu-Goto action 
reads 



<Sng — — 



2tt 



dv du£ , 



(8) 



where 



£ = y/1 - SIR 2 (QR' 2 + 2gj + / gg + ^,2) 



The equations of motion for <^(u) and R(u) follow from 
extremizing the Nambu-Goto action Q. One constant 
of the motion can be obtained by noting that the action 
is independent of f(u). This implies that 



n 
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(10) 



is constant. The minus sign on the right-hand side of 
( 10 1 is there in order to make II positive for positive 



f2. Eq. (10 1 can be solved for ip'(u) in terms of R(u) 
and R'(u). There are two solutions for f'{u): one in 
which ip'(u) is regular near the horizon at u = 1 and 
one in which tp'(u) diverges near the horizon. In infalling 
Eddington-Finkelstein coordinates the regular solution is 
aways the causal infalling solution and thus this is the one 
we choose. 

The equation of motion for R(u) is given by 



OR 



d <9£ 
dudR 1 



0. 



(11) 



Evaluating this expression and then eliminating (p deriva- 
tives via Eq. (10 1, we obtain the following equation of 
motion: 



R" 



R(u + 2RRf) (l + fRf 2 ) 

u(u 4 IL 2 - fR 2 ) 
u + ufR 12 -2(1-/) RR' (1 + Q 2 R 2 R' 2 ) 
uR(f-n 2 R 2 ) 



Eq. ( 12 ) is singular when 

u 4 IL 2 - fR 2 = or / - n 2 R 2 



(12) 
= 0. 

(13) 



As discussed in Refs. }45 , i;52j , reality of the Nambu-Goto 
action ([8| implies that both singularities must coincide at 
a single value of u, which we denote u = u c upon defining 



at which 
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R = R c = M-u c 



1/2 



(14) 



(15) 



Furthermore, goo = at the singular point u = u c . 
Consequently, there is a worldsheet horizon at u = u c , 
separating the upper part of the string u < u c which 
moves slower than the local velocity of light from the 
lower part of the string u c < u < Uh = 1 whose local 
velocity exceeds that of light. The string equation (12) 



becomes first order at the worldsheet horizon, 1 meaning 
that the equation of motion (12) itself determines R'(u) 
at u = m c , independent of any features of the solution 
away from (u c ,R c ). Consequently, II and 17 determine 
not only R(u c ) via Eqs. (14) and (15) but also R'(u c ), 



and then via the equation of motion ( 12 ) the entire solu- 
tion R(u) [52] , 

For a given II and £1 we determine R(u) numerically by 
solving Eq. (12) using pseudospectral methods [55]. We 
then solve Eq. (10) for <p(u), again using pseudospectral 
methods. 



IV. THE BULK TO BOUNDARY PROBLEM 

The presence of a string in the AdS-BH geometry, in 
our case the rotating string spiraling downward from the 
quark in circular motion at the boundary at u = 0, 
perturbs the geometry via Einstein's equations. In the 
N c — > oo limit, the (4+l)-dimensional gravitational con- 
stant becomes parametrically small and the presence of 
the string acts as a small perturbation on the AdS- 
BH geometry. To obtain leading order results in N c , 
we write the full (4+l)-dimensional metric as Gmn^= 



G 



(o) 

MN 



+ ^H MN , where G$ N is the AdS-BH metric |4|, 
and linearize Einstein's equations in the perturbation 
Hmn- This yields 



(16) 



where is a second order linear differential operator, 

k§ = 4n 2 /N 2 and 



t MN = __ 



-gg ab d a X M d b X N 5 i {r-r s ) (17) 



1 On can also show this directly by solving Eq. \l2\ via the Probe- 
nius method, i.e. doing a Laurent expansion about u = u c . 
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is the (4+l)-dimensional string stress tensor. The bound- 
ary value of the metric perturbation acts as a source for 
the J\f = 4 SYM stress tensor via the relation 1 60 



(18) 



H Llv =0 



where H flu (x) = lim M _>.o Hmn( x > u ) an d S grav is the grav- 
itational action. 

The (3+l)-dimensional stress tensor T^ v is symmetric, 
traceless and conserved and thus describes five indepen- 
dent degrees of freedom. These degrees of freedom can 
be conveniently packaged in five combinations of com- 
ponents of T M ", each of which transforms with definite 
hclicity under spatial rotations. There are two helicity 2 
components, two helicity 1 components and one helicity 
component of . Similarly, Hmn is a spin 2 field in 
(4+1) dimensions and thus contains 5 independent gauge 
invariant helicity degrees of freedom, gauge invariant in 
the sense that they are invariant under infinitesimal dif- 
fcomorphisms 



H 



MN 



Hmn — Dm£,n — -DwC 



AI- 



(19) 



where Dm is the covariant derivative with respect to the 
AdS-BH geometry and £m is an infinitesimal vector field. 
Each gauge invariant helicity degree of freedom in Hmn 
determines the corresponding helicity component of T^ v 
[551 [35]. Moreover, by rotational and gauge invariance, 
each gauge invariant helicity degree of freedom in Hmn 
satisfies a decoupled equation of motion. As we are only 
interested in the M — 4 SYM energy density, we focus 
below on helicity gauge invariants. 

To determine a helicity gauge invariant and its equa- 
tion of motion, it is convenient to decompose Hmn in 
terms of a complete set of functions \q ( x ) °f t ne 4<i co- 
ordinates x^. We choose 



Xq {x) = e-™ v i> q (x), 



(20) 



where tp q are eigenfunctions of the spatial Laplacian 
— V 2 ip q = q 2 ip q - In what follows, we shall initially not 
choose a specific basis of eigenfunctions ip q . We define 



MN 



J d 4 xx* q H M N, H qq = --^ J d 4 xx* q ViVjHij, 

(21) 

H q5 = — / d 4 XX* Q ^iH l5 , H 0q = ■ / d^XXn^iHoi, 

iq J q iq J q 

(22) 



and define similar expressions for tMN with the replace- 
ments H M n —> 1-mn an d Hmn — > t M N- (If the ip q were 
taken to be plane waves then Hmn would be the Fourier 
transform of Hmn-) Note that the q subscripts are la- 
bels, not values of the indices M or N, while the i and 
j are spatial indices corresponding to values of M or N 
given by 1, 2 or 3. When i or j is repeated, this indicates 



summation. With these definitions in hand, the helicity 
field that we shall use is given by 



Aq 2 



Hn 



MqfH 



(2uq 2 -3f) H qq 



(f-2uq 2 ) Hi 



<?5 



Hn 



(t 



ft 05 ) , 
(23) 



As can be easily verified, Z is invariant under the in- 
finitesimal diffcomorphisms (19 1. From the linearized 



Einstein equations (|16|) it is straightforward but tedious 

(24) 



to show that Z satisfies 

CZ = k 2 5 S, 
where the linear operator C is given by 



du 2 



q 2 (u 4 - 5) + 6u 2 (5u 4 - 9) 



- 2iu> 



u (q 2 + 6u 2 ) 

q 2 (5u 4 -5mw+9) - I8u 2 (3u 4 +3mw-7) - q 4 
u 2 (q 2 + 6u 2 ) 



d 
du 

2 



(25) 



and the source S is given by 
16(g 2 +18u 2 ) 

S = 8d u too 7 9 ; „ 9 x ^00 



4u (q 2 +6u 2 



u (q 2 +6u 2 



8iq(t 0q - ft q5 ) - 



8uq 2 f 
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3 

(2uq 2 - 



(tii ^t qq ) 



For strings which end at the boundary u 
source has the expansion 



3iw) t 05 . 

(26) 
0. the 



S = S {0) +S (1) u + O(u 2 ), 



(27) 



where 



'(0) 



-8 lim u 2 d u 

u-v0 



too 

1 1 



-8 lim d u to 



tt->0 



16 



2tt yr 



x* q (Ro) , 



(28) 



a result that follows from Eq. ( 26 ) upon noting that t, 



<<m> 



and indeed all the components of t in Eq. ( 26 1 are pro- 



portional to u for small u. Solving Eq. (24) with a series 



expansion near u = and demanding lim M _>o Hmn = 0, 
we find 



Z = Z(2)U 



(29) 



To lowest nontrivial order in u, Eq. (24 1 becomes just 
St, 



Z 



(2) 



'(0): 



but we shall also need because the 
change in the N = 4 SYM energy density due to the 
presence of the moving quark (i.e. the total energy den- 
sity minus that of the undisturbed plasma) is given by 



AS = - 



1 



(Z (3) +iujZ {2) ) . 



(30) 



(> 



This is as far as we can proceed without specifying a 
choice of the basis functions t\> q . 

We choose to compute the boundary theory energy 
density using the basis functions 



(31) 



where ji are spherical Bessel functions and Yi m are spher- 
ical harmonics. Our calculation is much simpler than 
the calculation for a quark moving along an arbitrary 
trajectory because, as in Ref. [45] . the quark has been 
moving in a circle at a constant angular velocity f2 for 
all time, meaning that everything in the problem (the 
4+1 dimensional string stress tensor and the 3+1 dimen- 
sional boundary theory energy density) is rotating at this 
constant angular velocity and depends on the azimuthal 
angle <j) and the time v only in the combination <p — Civ. 
(At the u — boundary, Eddington-Finkelstein time v 
is just boundary theory time t.) It follows that the en- 
ergy density only has support when u> = rn.fi. Hence, the 
Fourier transform in time is a discrete Fourier transform. 
Furthermore, specifying to specifies lo. 

For each {q, £, to} we solve Eq. ( 24 ) using pseudospec- 



tral methods in which the u-dependence of all functions 
is decomposed in terms of Chebyshev polynomials |59j . 
The requisite two boundary conditions are specified at 
U = and u = 1, as we now describe. At u = 0, there 
are two linearly independent solutions. The one we want 
is regular at u = and has the expansion ( 29 ) , while 



the one that must be taken to vanish is proportional to 
it 3 log w at small u. Imposing a boundary condition that 
selects the regular solution is made easier by first writing 



Z = u 2 S t 



(0) 



u A X 



(32) 



and turning the differential equation for Z into a differ- 
ential equation for X. In this differential equation, X" 
(by ' we mean d u ) arises only in the term ufX", mean- 
ing that at ii — O the differential equation for X involves 
only X and X', as long as X" is finite. So, the differen- 
tial equation itself specifies its own boundary condition 
(as a relation between X and X') at u = 0. Satisfying 
this boundary condition automatically yields a regular 
solution at u = 0, since if X is finite the boundary con- 
dition makes X' finite, and if X' is finite the boundary 
condition makes X" finite, and so on. (The undesirable 
solution has X, X' and X" all divergent at u = 0.) At 
the AdS-BH horizon u = 1, there are again two linearly 
independent solutions and we must specify a boundary 
condition that selects only the infalling mode. Because 
we are using ingoing Eddington-Finkelstein coordinates, 
this boundary condition is equivalent to requiring that Z 
(or X) and all of its derivatives are finite at u = 1. Again, 
as long as X" is finite at u = 1 the term ufX" vanishes 
there, this time because / = at u = 1, and the dif- 
ferential equation for X turns into a boundary condition 
relating X and X' which, when satisfied, yields a regu- 
lar solution at u = 1. (In ingoing Eddington-Finkelstein 
coordinates, the outgoing mode has a divergent phase at 
u = 1.) 



Once the solution Z(q, i, to, u) is determined, we ex- 
tract Z( 3 ) (which is just X evaluated at u = 0) and 
compute the Fourier components of the energy density 
A£(q,£,m) via (30). The real-space energy density is 
then computing by evaluating 

,( 33 ) 

numerically. This is the change in the energy density at 
time t = relative to the energy density of the unper- 
turbed plasma due to the presence of the rotating quark. 
The time-dependent energy density A£(r, 6, </>, i) is ob- 
tained simply by replacing <f> by <j) — fit. In order to ob- 
tain the results that we illustrate in the next Section, we 
typically used between 10 3 and 10 4 values of q between 
and 200 irT and at least 40 values of £, with smooth 



window functions cutting the Fourier transform ( 33 ) off 
at large q and I. 

One further complication remains, before we turn to 
our results. The differential operator C and the source S 



in ( 25 ) and ( 26 ) are badly behaved at u = for small q. 



There is no problem of principle, but solving the equa- 
tion (24) numerically becomes intractable. Because the 



energy density has support only where g > uj = mfl, this 
difficulty only arises for the to = modes. We therefore 
use the procedure that we have described only for the 
to +~ components of the energy density A£(g, £, to). We 
calculate the m = components using a different gauge 
invariant helicity field Z that we define and describe 
in Appendix A. 



V. RESULTS 

Fig. [l] shows three different plots of r 2 A£ / P for quarks 
in circular motion with each of three different velocities: 
/3 = 0.15, $ = 0.3 and /3 = 0.5. Here, P = dE/dt is the 
energy lost by the circulating quark (and hence dumped 
into the plasma) per unit time. The radius of the quark's 
trajectory in all plots is Rq = 1 and the temperature of 
the plasma is given by 7rT = 0.15/i?o- This means that 
S defined in is given by 1.0, 4.6 and 17.1 in the left, 
middle and right columns respectively. P is given by 
([2]) when S ^ 1 or by Q when S « 1, P has been 
calculated in Ref. [52 for any S and is related to the 
constant n defined in (fl0| by 



P 



2tt 



-mi 



v/Xn/3 



(34) 



In the left, middle and right columns of Fig. [TJ n is given 
by 5.2 x 10- 3 /i? , 3.0 x lO" 2 /^ and 0.22/i? respec- 
tively. At the time shown, the quark is located at x = Rq, 
y = and the quark is rotating counter-clockwise in the 
plane 2 = 0. The three plots in the top row are cut- 
away plots with the cutaways coinciding with the planes 
z = 0, <f) = and <\> — 7tt/5. The three plots in the 
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FIG. 1: Plots illustrating the energy density of strongly coupled M = 4 SYM plasma in which a test quark is rotating on a 
circle with radius Rq with angular velocity Q = /3/Ro for /3 = 0.15 (left column), /3 = 0.3 (middle column) and j3 = 0.5 (right 
column), corresponding to H = 1.0, 4.6 and 17.1. In all plots, the temperature of the plasma is given by irT = 0.15/i?o and the 
units are chosen such that Ro = 1. Top: cutaway plots of r 2 A£/P where P is the power radiated by the quark. The cutaways 
coincide with the planes 8 = 7r/2, <f> — and 4> — 7ty/5. Middle: plots of r 2 A£/P on the equatorial plane 9 = 7r/2 (i.e. z = 0). 
Bottom: blue curves are plots of r 2 A£/P at 6 — ty/2 and cf> = it/2. The quark's trajectory lies in the equatorial plane 6 — ty/2 
and the quark is rotating counter-clockwise. The red dashed curves in the bottom plots show r 2 £/P for the strongly coupled 
synchrotron radiation emitted by a quark in circular motion in vacuum [45], pulses of radiation that propagate outward to 
r — > oo at the speed of light without spreading. 



middle row show the energy density on the plane z — 
and the bottom three plots give the energy density at 
z = 0, = 7r/2, namely a slice through the middle-row 
plot along one radial line. For reference, the red dashed 
curves in these bottom plots show r 2 £ for the strongly 
coupled synchrotron radiation that a quark moving along 
the same circular trajectory would emit in vacuum |45j . 
In each of the bottom plots, we use the same P to nor- 
malize the red curve as for the blue curve. All nine panels 
in Fig. [I] show the energy density at one instant of time, 
but the time-dependence is easily restored by replacing 
the azimuthal angle <fi by <fi — f2i, where — ft/Ro is 
the angular velocity. As a function of increasing time, 
the entire patterns in the upper and middle rows rotate 
with angular velocity f2, as the spirals of radiation move 
outwards. As a function of increasing time, the patterns 
in the lower rows move outwards, repeating themselves 
after a time 2ir/il. 

As is evident from Fig. [lj as the quark accelerates 
along its circular trajectory, energy is radiated outwards 
in a spiral pattern which is attenuated as the radia- 
tion propagates outwards through the plasma to increas- 
ing r. However, the qualitative features of the spiral 
patterns differ greatly at the three different quark ve- 
locities that we have chosen. For j3 = 0.15 the spi- 
ral arms are very broad in r, as broad as their separa- 
tion, and the spiral pattern propagates outwards at the 
speed of sound, while being attenuated with increasing r. 
(Second-order hydrodynamics for a conformal fluid with 
a gravity dual like Af — 4 SYM theory predicts a sound 
velocity 1/V3 + 0.116 q 2 /{irT) 2 + . . . [BT] for sound waves 
with wave vector q. The sound waves in the left column 
of Fig. [TJhave q ~ 1.3 ttT and are moving outward with 
a velocity of 0.73. We shall return to the comparison 
to second-order hydrodynamics in Section VI.) The red 
curve in the lower-left panel shows the energy density 
of the synchrotron radiation that this quark would have 
emitted if it were in vacuum, and we see that there is no 
sign of this in our results. So, at this /?, corresponding 
to S = 1.0, the rotating quark is emitting sound waves. 

The results in the right column of Fig. [T] for f3 = 0.5, 
are strikingly different. The spiral arms are very nar- 
row in r, much narrower than their separation, and they 
propagate outwards at the speed of light, as can be seen 
immediately in the bottom-right panel by comparing our 
results, in blue, to the energy density of the synchrotron 
radiation that this quark would have emitted if it were 
in vacuum. We see that at this /3, corresponding to 
S = 17.1, the rotating quark is emitting strongly cou- 
pled synchrotron radiation, as in vacuum |45j . and we 
see that the radiation is being attenuated as it propagates 
outward in r, through the strongly coupled plasma. Re- 
markably, even as the outgoing pulses of energy are very 
significantly attenuated by the medium we see no sign of 
their broadening in either the 9 or the <j> or the r direc- 
tions. Looking at the vertical sections in the upper-right 
panel, we see that if anything the spread of the beam of 
radiation in is becoming less as it propagates and gets 




FIG. 2: Energy density r 2 £/P of the strongly coupled syn- 
chrotron radiation emitted by a quark in circular motion with 
Ro — 1 and j3 — 0.5, exactly as in the top-right panel of Fig. IT] 
but this time in vacuum, at T = 0, calculated as in Ref. [45]. 
The color-scale is the same as in the top-right panel of Fig.[T] 
and the red curve in the bottom-right panel of Fig. [l] is the 
profile along one radial line through this figure. 



attenuated. This conclusion is further strengthened by 
comparing the upper-right panel of Fig. [I] to the anal- 
ogous results for a quark in circular motion in vacuum 
shown in Fig. [2] It is certainly clear that the presence 
of the medium does not result in the spreading of energy 
away from the center of the beam at the equator out 
toward large polar angles. Just the opposite, in fact: at 
large polar angles the beam gets attenuated more rapidly 
than near 8 = tt/2. We shall return to this in Section VI. 
Broadening in either the <f> or the r directions would be 
manifest as widening of the pulses in the bottom-right 
panel, and this is certainly not seen. In fact, we have 
extended the plot in the bottom-right panel out to larger 
r, for several more turns of the spiral, and we continue 
to see rapid attenuation with no visible broadening. 

We turn our attention now to the center column of 
Fig. [T] Here, with a rotation velocity of j3 = 0.3 cor- 
responding to E = 4.6, we clearly see both synchrotron 
radiation and sound waves. The synchrotron radiation is 
most easily identified with reference to the results for a 
quark with this rotation velocity in vacuum, shown in the 
red curve in the bottom-center panel. In our results with 
T 0, we see the emission of a pulse of synchrotron radi- 
ation whose amplitude is very rapidly attenuated, much 
more rapidly than in the right column. In part guided by 
our inspection of the results at large H in the right col- 



= 0.65 



umn, we see that as the pulse of synchrotron radiation is 
attenuated, it too does not broaden. What we see here 
that is not so easily seen in the right column is that as the 
pulse of synchrotron radiation is attenuated it "sheds" a 
sound wave, leaving behind it a broad wave, reminiscent 
of the sound waves in the left column. Behind each pulse 
of synchrotron radiation we see the "compression half" 
of a sound wave, and behind that a deeper rarefaction, 
and then the next pulse of synchrotron radiation arrives. 
Once seen in the middle column, this phenomenon can 
perhaps also be discerned to a much lesser degree in the 
right column, with the each pulse of synchrotron radia- 
tion trailed first by a slightly yellow region of compres- 
sion and then by a more blue region of rarefaction. It is 
not really clear in the right column whether these can be 
called sound waves, both because of their smaller ampli- 
tude and because the next pulse of synchrotron radiation 
overwhelms them sooner than in the middle column. In 
the middle column, though, the interpretation is clear: 
the beam of synchrotron gluons is exciting sound waves 
in the plasma. In the right column, it is clear that any 
sound waves, if present, are smaller in amplitude than in 
the middle column even though the pulses of radiation 
carry more energy in the right column. 

Fig. [^demonstrates clearly that at small S the rotating 
quark emits only sound waves while at large S it emits 
strongly coupled synchrotron radiation as in vacuum, 
with that beam of gluons subsequently being quenched 
by the plasma. The crossover between these two regimes 
lies between S = 1 and S = 5, consistent with expecta- 
tions based upon the results of Ref. [52]. We have also 
confirmed that this crossover occurs in the same range 
of S for a quark in circular motion with RqttT twice as 
large as in Fig. [l] meaning that this range of S occurs at 
larger (3. For RqitT = 0.15 as in Fig. [TJ we have looked 
at a quark moving with ft = 0.10, meaning S = 0.45 
and confirmed that we see only sound waves, with longer 
wavelengths than in the left column of Fig. [I] and hence 
with a velocity that is closer to l/v3, the q — > ve- 
locity of sound. We have also done the calculation at 
P = 0.65, corresponding to 3 = 42.8, with results shown 
in Fig. [3] which we have extended out to larger r to show 
the beam of synchrotron radiation getting almost com- 
pletely attenuated. The comparison in the lower panel 
between the pulses of radiation propagating through the 
plasma (blue curve) and those propagating in vacuum 
(red dashed curve) makes it clear that even as the beam 
is being almost completely attenuated by the plasma, it 
propagates at the speed of light and, as far as we can see, 
it does not broaden. 

Extending our calculations with RqttT = 0.15 to larger 
7 is possible, but the numerics rapidly become more dif- 
ficult as the radial width of the pulses narrows like I/7 3 , 
rapidly increasing the required dynamical range in mo- 
mentum space and rapidly making the Fourier transform 
back to position space more costly. If it were impor- 
tant to pursue this, however, it could certainly be done. 
But, we shall see in Section VI that we have an analytic 
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FIG. 3: Energy density r 2 A£/P of strongly coupled Af = 4 
SYM plasma in which a quark is rotating on a circle with 
RqttT = 0.15 and with velocity f3 = 0.65, meaning H = 42.8. 
Plots as in one column of Fig. [I] extended to r = 120 to show 
the quenching of the beam of strongly coupled radiation. 



understanding of the qualitative features in our results, 
and that based upon this analytic understanding we do 
not expect any qualitatively new behavior at larger 7 at 
the same Rq. However, it would be very interesting to 
explore the regime in which RqttT 3> 1, meaning that 
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Fl/irT <C 1, and yet 7 is so large that S > 1. In this 
regime, we expect narrow pulses of radiation that are 
much more widely separated in the radial direction than 
those we have analyzed, say in the right column of Fig. [T] 
This would allow us to learn more about the response of 
the strongly coupled plasma to a beam of synchrotron 
radiation since we would be able to watch the plasma 
for much longer after one pulse of radiation passed be- 
fore the next pulse arrived. Unfortunately, there is at 
present a serious obstacle to doing such a calculation, 
more serious than can be overcome simply by increas- 
ing the dynamical range of the calculation. We saw in 
Section IV that the helicity gauge invariant quantity 
Z that we defined in Section IV has gravitational equa- 
tions of motion that are badly behaved for q <C 7rT, but 
that because the energy density has support only where 
q > lo — mO, with m an integer, this problem only arose 
for tyl — modes — which we were able to analyze us- 
ing a different gauge invariant described in Appendix A. 
However, if f2 -C irT then there are many modes with 
u) 7^ and uj <C ttT which are important. We have not 
found a gauge invariant that yields a tractable numerical 
analysis for such modes — the analysis of Appendix A 
only works for u) = 0. 

There are in fact several (related) obstacles to using 
our calculation to provide definitive answers to the ques- 
tion of where the energy that is initially in the gluon 
beam goes as the gluon beam gets attenuated. The first 
we have discussed above: we cannot watch the plasma 
behind one of the pulses of radiation very long before 
the next pulse comes along and obliterates whatever the 
previous pulse has left behind. However, even if a reso- 
lution to the technical obstacle that currently precludes 
addressing this were found, a further obstacle remains. 
We are analyzing a scenario in which the quark has been 
moving in a circle for an infinitely long time meaning that 
a steady-state in which the energy density at any posi- 
tion is a periodic function of time has been achieved. We 
see in Figs. [T] and [3] that the energy density in the beam 
falls off faster than 1/r 2 at large r. So, the natural first 
expectation is that the beam heats the plasma up in the 
range of r over which it gets attenuated — perhaps it first 
makes sound waves, but ultimately these too will damp, 
leaving just a heated region of plasma. This expectation 
cannot be correct in a steady-state calculation like the 
one we have done, since a continual heating up of some 
region of space blatantly contradicts the steady-state as- 
sumption. So, what actually happens to the energy in 
our calculation? We have checked that at sufficiently 
large r the energy density A£ is zero. This means that 
at sufficiently large r, there is an outward flux of energy 
whose magnitude, averaged over angles, is P/(4irr 2 ) with 
P the energy lost by the rotating quark per unit time. 
This energy flux corresponds to a collective outward flow 
of the plasma with a velocity, averaged over angles, given 
by 



where we have used the fact that the sum of the en- 
ergy density and pressure of the plasma in equilibrium 
is £ + p = 7r 2 A^T 4 /2. Since we are working in the 
large- N c limit, the velocity v p i aS ma is infinitesimal. So, in 
our steady-state calculation, the energy from the gluon 
beam ulltimately finds its way into an infinite wavelength 
mode with infinitesimal amplitude. A mode like this can 
equally well be thought of as a sound wave with infinite 
wavelength and infinitesimal amplitude (i.e. infinitesimal 
longitudinal velocity) or as a diffusive mode with infinite 
wavelength. 2 This is the only possible answer to the ques- 
tion of where the energy from the gluon beam ultimately 
ends up in a steady-state calculation like the one that 
we have done. In a sense, this energy flux corresponding 
to an infinitesimal- velocity outward flow of the plasma is 
the closest that a steady-state calculation can come to 
describing the heating up of a region of the plasma — 
which cannot happen in steady-state. 



VI. DISCUSSION 

We turn now to a discussion of our results. Much can 
be learned about the qualitative features of the results 
illustrated in Fig. [I] by studying the quasinormal modes 
of the AdS-BH spacetime that provides the dual gravita- 
tional description of the physics. In the dual gravitational 
picture, the moving string excites a full spectrum of grav- 
itational quasinormal modes, which propagate outwards 
and eventually get absorbed by the black hole. The prop- 
agation and absorption of these quasinormal modes man- 
ifests itself on the boundary as the propagation and at- 
tenuation of the spirals of energy density shown in Fig.[T] 
The dispersion relations uj(q) of the helicity quasinor- 
mal modes are given by solutions to [52] 



det £(uj,q) = , 



(36) 



where the linear operator £(uj,q) is defined in Eq. (25) 



1 



fpla 



47rr 2 (£ + p) 2N? (ttT) 2 (rirT) 2 



(35) 



and where the determinant is to be evaluated in the space 
of functions satisfying the appropriate boundary condi- 
tions that we have described earlier. Upon introducing a 
complete basis of functions of u that satisfy the boundary 
conditions and then truncating that basis, C becomes a 
matrix and Eq. (36) can be solved for uj{q) numerically. 

Fig. [4] shows a plot of the quasinormal mode disper- 
sion relation for the lowest quasinormal mode (i.e. the 
one with the smallest imaginary part). For q <C ttT the 
dispersion relation has the asymptotic form expected for 



2 In a relativistic plasma, an infinitesimal increase in the tempera- 
ture in some region must correspond to an infinitesimal increase 
in the pressure in that region, meaning that it corresponds to 
sound waves. So, in this case this mode is better thought of 
as an infinite wavelength sound wave rather than as a diffusive 
mode. 
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FIG. 4: A plot of the real and imaginary parts of the dis- 
persion relation of the lowest quasinormal mode. We plot 
Reuj/q and Imu / '{irT) since these ratios are both of order 1. 
For q <C nT the dispersion relation is that of sound wave s 
whose dispersion relation is given up to order q 3 by Eq. (37 1, 
plotted as dashed lines in the figure. For q ;§> nT the disper- 
sion relation is that of waves which propagate at the speed 
of light. The large-g asymptotic expression (39 1 that we have 



obtained by fitting the results in this figure is plotted as the 
dotted lines. 



the hydrodynamics of any conformal fluid [61) 



u s (q) = c s q - 



iTq 2 + - 



2' 9 



+ 0(q 4 ) , (37) 



where in Af = 4 SYM theory, with its classical gravity 
dual, all the constants are known analytically: the low- 
q speed of sound is c s = I/a/3, the sound attenuation 
constant T is given by ttTT = 1/6, and the relaxation 
time Tn is given by 7rTr n = (2 — log2)/2. These modes 
represent propagating sound waves which attenuate over 
a timescale 



^damping 



1 



(38) 



The dispersion relation (37) is plotted in Fig. |4j it de- 
scribes the full dispersion relation very well for q < 2wT. 
This supports our observation that the waves in the left 
column of Fig. [I] are sound waves. Since these waves are 
not monochromatic (and since in the dual gravitational 
description they are not described solely by the lowest 
quasinormal mode) they cannot be compared quantita- 



tively to (37), but their velocity is as (37) predicts for 
which is comparable to the q ~ 1.35 nT 



q ~ 1.2 nT, which is comparable to t: 
obtained from their peak-to-peak wavelength. Using 
q ~ 1.2 7rT in ( 38 ) predicts a sound attenuation timescale 
(Ai) soun( j ~ 4^5/ (7rT) ~ 30i?o, which is comparable to 
but a little shorter than the exponential decay time for 
the amplitude of the waves in the left column of Fig. [T] 
which is closer to 40 Rq. So, although a quantitative com- 
parison is not possible, it does seem that the low-q regime 



of the dispersion relation in Fig. [4] that describes sound 
waves does a reasonable job of capturing the qualitative 
features of the waves seen in the left column of Fig. [T] 

The dispersion relations of the higher quasinormal 
modes (those with more negative imaginary parts) can 
also be determined by solving Eq. ( 36 ) numerically. At 
q <C nT they approach the asymptotic form lj = (a — 
ib)nT where a and 6 are mode-dependent 0(1) con- 
stants, with values that are larger and larger for higher 
and higher modes. (For the lowest quasinormal mode, 
5 = 6 = 0.) At low q, disturbances of the plasma 
described by higher quasinormal modes attenuate on a 
timescale of order 1/(6 nT) that is much shorter than 
that for the sound waves described by the lowest quasi- 



normal mode, namely (38) 



Let us turn now to q ^> nT. By fitting to our results for 
the dispersion relation for the lowest quasinormal mode, 
we find that in this regime the dispersion relation takes 
the asymptotic form 



' 7rT\ 

^rad = Q + KT(a - lb) ( — J 



(39) 



as argued for on general grounds in Ref. [S3], where we 
find a « 0.58 and 6 sa 1.022. At q > 7rT the dispersion 
relations of all quasinormal modes approach the asymp- 
totic form (39), with a and 6 mode-dependent 0(1) con- 
stants, again with values that are larger and larger for 
higher and higher modes. Therefore, generically the high 
q modes propagate at close to the speed of light and at- 
tenuate over a time-scale 



t 



damping 
fad 



l 



irTb \ttTJ 



1/3 



(40) 



where we shall use the value 6 « 1.022 from the lowest 
quasinormal mode in making estimates, keeping in mind 
that if the contribution of higher quasinormal modes were 
important this would increase the effective 6 somewhat. 
The fact that the pulses of energy in Fig. [T] are far from 
being monochromatic waves introduces a larger uncer- 
tainty into our discussion than does not knowing how 
much the higher quasinormal modes contribute. 



We have plotted the large-g asymptotic expression ( 39 ) 
for the dispersion relation for the lowest quasinormal 
mode in Fig. |1J and we see that it describes the full re- 
sult very well for q > 20 7rT, and has the right shape 
at a qualitative level down to about q ~ 5 ttT. This is 
consistent with our observation that the narrow pulses 
of synchrotron radiation in the middle column, where 
the pulses have a full width at half maximum (FWHM) 
~ 2.5 i?o corresponding very roughly to q ~ 6 7rT, and 
the right column, where the pulses have a FWHM ~ i?o 
corresponding very roughly to q ~ 15 7rT, propagate out- 
wards at the speed of light. Converting the widths of 
these pulses into estimates of q is very rough because the 
pulses are neither sinusoidal nor Gaussian. If we never- 
theless try substituting q ~ 15 7rT into (40) we find that 
it predicts t^™ pins ~ 16i?o, which is roughly half the 
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exponential decay time for the amplitude of the waves 
in the lower-right panel of Fig. [TJ Again, quantitative 
comparison is not possible, but inferences drawn from 
the large-g dispersion relation for the lowest quasinormal 
mode ( 39 1 is at least in the right ballpark. 



the strongly coupled plasma itself and have nothing to do 
with the source on the right-hand side of ( 24 1 , namely the 



We can also compare t^™ pins for the radiation emitted 
by the quark moving with j3 = 0.5 in the right column 
of Fig. [T] with that emitted by the quark moving with 
/3 = 0.65 in Fig. [3] Since the pulses of synchtron radia- 
tion narrow proportional to I/7 3 as (3 increases, Eq. (40) 



predicts that t^™ pmg should increase proportional to 7, 
namely by about 14% in going from f3 = 0.5 to (3 = 0.65. 
It is hard to define and extract a damping time from 
the figures with sufficient accuracy to test this prediction 
quantitatively but, again, it is in the right ballpark. 



The qualitative prediction from ( 40 ) that narrower 
pulses, with higher q, can penetrate farther into the 
strongly coupled quark-gluon plasma is apparent in com- 
paring the ^-dependence of the results in the top-right 
panel of Fig. [I] with the T = results in Fig. [5j We 
see that in vacuum the pulses of synchrotron radiation 
get broader and broader as you decrease 9 from 7r/2, 
going from the equator toward the poles. And, indeed, 
just as (40) predicts we see in Fig. [I] that in the quark- 
gluon plasma the radiation heading in polar directions is 
quenched much more quickly than that radiated in nar- 
rower pulses near the equator. 

We can also use the quasinormal mode dispersion re- 
lation to understand why the pulses do not broaden sig- 
nificantly in the radial direction as they propagate. The 
increase in the width of a pulse as it propagates for a 
time t is ~ t AqcPui/dq 2 , where Aq is the width of the 
pulse in q-space. Taking Aq ~ q and using the large-g 
dispersion relation ( 39 ) , we find that after the radiation 



damping time given by ( 40 1 the pulse should have broad- 
ened by ~ 4a/(9bq). If the pulse had a Gaussian profile, 
this would correspond to broadening by about 10% of 
the original FWHM of the pulse. So, the quasinormal 
mode dispersion relation predicts that by the time the 
pulses have been significantly attenuated, they should 
have broadened by an amount that is parametrically of 
order their initial width, but smaller by a significant nu- 
merical factor. It is therefore not surprising that we see 
no significant broadening in Fig. [T] 

By this point we have understood many of the most 
interesting features of our results qualitatively, and even 
scmi-quantitativcly, by analyzing the quasinormal mode 
dispersion relations. This gives us confidence that no 
new qualitative phenomena emerge for narrower pulses 
(higher g; e.g. from a rotating quark with larger 7) than 
we have been able to study, since it is clear that our 
results at (3 — 0.5 and (3 = 0.65 are already exploring 
the high-g regime of the dispersion relation in Fig. |4j 
where the asymptotic expression (39) is a good guide. It 



is also important to stress that the quasinormal mode 
frequencies are determined entirely by £, from the left- 



rotating quark. Given that we have been able to use the 
quasinormal mode dispersion relations so successfully to 
understand the propagation, rate of attenuation and lack 
of broadening of a beam of gluons, we are confident that 
these phenomena are are independent of how the beam 
of gluons is created. 

Finally, we can use our understanding of the quasinor- 
mal mode dispersion relations to clarify the distinction 
between the middle column of Fig. [T] — in which we see 
the pulse of radiation shedding a sound wave — and the 
cases where we don't see this (like the right column of 
Fig. [l] and Fig. [3|) As we described at the end of Section 
V, the fact that we have done a steady-state calculation 
makes it difficult to use our results directly to draw con- 
clusions about what our pulse of radiation leaves behind. 
However, we can use our understanding of the quasinor- 
mal modes to answer this question at least qualitatively. 
Suppose that we could move the quark on some trajec- 
tory such that it radiates one pulse of synchrotron radia- 
tion (i.e. one turn of the spiral) and then no more; what 
would happen to this single pulse of strongly coupled ra- 
diation as it propagates outward through the strongly 
coupled plasma? The dual gravitational description of 
this radiation would be governed by ( 24 ) , with the same 



C and hence the same quasinormal modes as in our anal- 
ysis but with a different source S, localized along the 
world sheet of the string hanging down from the quark 
that excited the single pulse of radiation. As long as we 
look only at distances greater than of order 1/(ttT) away 
from the location of the string, the disturbance of the 
plasma must be described by a pulse of short wavelength 
radiation with the dispersion relation ( 39 ) that moves at 



the speed of light and is attenuated on timescales (40) 



as well as long wavelength sound waves with the disper- 
sion relation (37) that propagate outward at the speed 
of sound and are attenuated on timescales (38). Since 



these sound waves move more slowly, the pulse of radi- 
ation leaves them behind — shedding them as we see in 
the middle column of Fig. [T] The same would happen for 
shorter wavelength pulses as in the right column of Fig. [T] 
or as in Fig. [3J but in these cases in our steady-state cal- 
culation the next pulse of synchrotron radiation arrives 
before we can see the sound waves being left behind. In 
the hypothetical case of a single short wavelength pulse, 
the short wavelength pulse itself will get far ahead of the 
sound waves it has left behind as it is attenuated only 
on the long timescale ( 40 1 . By the time the short wave- 



hand side of (24). This means they reflect properties of 



length pulse has damped away, the sound waves that it 
shed will be far behind and will have become sound waves 
with very small q, meaning long wavelengths and small 
amplitudes, since (38) tells us that these are the sound 
modes that last the longest. These are represented in 
the steady-state calculation by the outward-going energy 
flux at infinite wavelength and infinitesimal amplitude 
that we described at the end of Section V. We now see 
that the distinction between the middle column of Fig. [I] 
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on the one hand and the right column of Fig. [T] and Fig. [3] 
on the other hand is that in the former case the pulse of 
radiation is never well-separated from the sound waves 
that it leaves behind — if a significant pulse of radiation 
remains, the sound waves are not far behind it and have 
themselves not yet been attenuated — while in the latter 
case the pulse of radiation gets far ahead of the sound 
waves and the sound waves "thermalize" (which means 
increase in wavelength and decrease in amplitude) before 
the pulse of radiation has been attenuated. 



VII. FROM QUENCHING A BEAM OF 
STRONGLY COUPLED GLUONS TO JET 
QUENCHING 

There are many qualitative similarities between the 
quenching of the beam of strongly coupled synchrotron 
radiation in the strongly coupled J\f = 4 SYM plasma 
that we have studied and the quenching of jets in heavy 
ion collisions at the LHC and RHIC: 

• As our beam of gluons propagates through the 
plasma, losing a significant fraction of its energy, 
it does not spread in angle. Jets in heavy ion colli- 
sions at the LHC lose a significant fraction of their 
energy but, within current experimental errors, do 
not broaden in angle and do not get deflected in 
their direction. 

• As our beam of gluons propagates through the 
plasma, and is significantly attenuated, it does not 
spread in the direction along which it propagates. 
In other words, the momenta of the gluons making 
up the beam do not decrease even as the beam loses 
a significant fraction of its energy. Similarly, jets in 
heavy ion collisions at the LHC have fragmenta- 
tion functions (i.e. distributions of the momenta of 
the particles making up the jet) that are unmod- 
ified by propagation through the strongly coupled 
plasma produced in the collisions except by virtue 
of the overall reduction in the energy of the jet. 

• In the case of a beam made up of gluons whose 
wavelength is not too short, as in the middle col- 
umn of Fig. [TJ it is possible for a significantly atten- 
uated pulse of radiation to be followed close on its 
heels by a significant sound wave, trailing behind it. 
In contrast, once a beam made up of shorter wave- 
length gluons, as in the right column of Fig.[T]or in 
Fig. [3j is significantly attenuated, the sound waves 
that it left behind are far behind it and have them- 
selves thermalized. In our steady-state calculation, 
they have ended up in the infinite wavelength, in- 
finitesimal velocity, mode that we described at the 
end of Section V. The behavior of these shorter 
wavelength pulses, losing their energy to a mode 
in which the collective velocity of the fluid is in- 
finitesimal, suggests the observation that jets at the 



LHC lose their energy to soft particles at all angles 
relative to the jet direction. The behavior of the 
longer wavelength pulses suggests that there may 
be a regime of jet energies, analogous to the beam 
of radiation in the middle column in Fig.[l] in which 
at a time when the jet itself has been attenuated 
significantly it is followed by a significant pulse of 
sound waves moving in the same direction as the 
jet. If the jets accessible in RHIC collisions, with 
energies in the 20-40 GeV range, are in this regime, 
this would indicate that they should lose their en- 
ergy to soft particles that are correlated in angle 
with the jet direction, consistent with preliminary 
indications from RHIC data. 

Comparisons along these lines will never be more than 
qualitative, since the beam of strongly coupled radia- 
tion that we have analyzed is not a jet. However, these 
multiple qualitative resonances between jet quenching in 
heavy ion collisions and the quenching of the beam of 
strongly coupled radiation that we have analyzed support 
the prospect that jet quenching is a strongly coupled phe- 
nomenon, even for the few-hundred-GeV jets produced at 
the LHC. If this is so, what is to be done next? Further 
directions include: 

• Finding a helicity gauge invariant that makes 
the regime with R ttT ^S> 1 and 7 large enough 
that S 3> 1 accessible, since this would give much 
more time for the region of plasma through which a 
pulse of radiation has passed to be analyzed before 
the next pulse hits and thus would permit a more 
definitive analysis of where the energy lost by the 
quenched beam ends up. 

• Finding ways to make a beam of strongly coupled 
radiation other than via synchrotron radiation, and 
in particular finding a way to make such a beam 
pointing in a pair of fixed directions rather than 
rotating. For example, could multiple quarks mov- 
ing in concert be engineered so as to behave like a 
phased array of antennas, generating back-to-back 
beams of strongly coupled radiation? 

• Reanalyzing the stopping of a light quark or gluon, 
as in Refs. [40l - l42| . should give another means of ac- 
cessing many of the questions we have addressed. 
And, these analyses find a stopping distance for 
energetic light quarks that is proportional to E 1 / 3, , 
with E the energy of the quark. Our result for 
the distance scale over which our beam of gluons is 
quenched, namely t^™ pm& ~ q 1 / 3 j (7rT) 4 / 3 where 
q is the typical wave vector of the gluons in the 
beam, has the same parametric dependence as the 
light-quark stopping distance, adding considerable 
robustness to both results. There have been some 
first steps taken to compare this relationship be- 
tween quark or gluon energy or wave vector and 
stopping or quenching distance with heavy ion col- 
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lision data | 
these efforts. 



Our results further motivate 



• It is also interesting to note that in the calculation 
of Ref. [32] a very high energy quark that loses, 
say, half of its energy leaves most of that energy 
far behind it while a lower energy quark that loses 
the same amount of energy and comes almost to 
rest is never well-separated from the energy it has 
lost. Although described in quite different terms, 
this is reminiscent of the distinction between the 
middle and right panels of Fig. [T] and perhaps of 
the possible distinction between jet quenching at 
the LHC and at RHIC. This distinction may also 
be seen in the analysis of the energy loss of a heavy 
quark being dragged through the strongly coupled 
plasma [37J UHL m which the world sheet of the 
string trailing behind the dragging quark features a 
world-sheet horizon [2"9"rf3"T] , located near the space- 
time horizon for a low velocity quark and located 
closer and closer to the boundary for quarks mov- 
ing more and more relativistically. One possible in- 
terpretation of the world-sheet horizon is that the 
portion of the string between it and the boundary 
describes (ultraviolet) modes in the gauge theory 
that propagate along with the heavy quark while 
the portion of the string between the world-sheet 
horizon and the spacetime horizon describe the dis- 
turbance of the plasma that the quark leaves be- 
hind [66-69]. This interpretation predicts that a 
higher energy heavy quark will be accompanied 
only by shorter wavelength modes of the gauge the- 
ory while a lower energy heavy quark will be accom- 
panied by softer gauge theory modes. Although de- 
scribed in terms that are different yet again, this is 
again reminiscent of the possibility that only lower 
energy quenched jets will be accompanied by soft 
particles correlated in angle with the jet direction. 
All these on-the-surface quite distinct approaches 
to jet quenching point in the same direction, sug- 
gesting that the distinction manifest in Fig. [T] be- 
tween an attenuated high energy jet that has left its 
sound waves far behind and an attenuated lower en- 
ergy jet that has a pulse of sound close on its heels 
may be generic to jet quenching in strongly coupled 
plasma. The connections between these different 
approaches must be explored and developed. 

• Our results also motivate further effort to deter- 
mine whether the soft particles in a heavy ion colli- 
sion corresponding to the energy lost by a quenched 
jet really are more correlated with the jet direction 
at lower jet energies. If so, this will point to lower 
jet energies as the place to look for the hydrody- 
namic response of the plasma. 

• On a more theoretical note, it would be very in- 
teresting to repeat our analyses in a nonconformal 
strongly coupled gauge theory plasma with a dual 
gravitational description. 
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Appendix A: Zero Mode Gauge Invariant 

The helicity gauge invariant Z that we defined in 
Section IV is inconvenient at m = lj = because the 
equations for it have apparent divergences that preclude 
their numerical solution. In the static case with w = 0a 
different helicity gauge invariant 

Z =q 2 n 00 + y (2 - /) (Hu - H qq ) (Al) 

proves convenient instead. 3 As one can easily verify, Z 
is invariant under time-independent infinitesimal diffeo- 
morphisms and satisfies the equation of motion 

CoZ = kIS , (A2) 

where the linear operator C is given by 

d 2 9-8u 4 + 7u 8 d 



^ du 2 u(3 — it 4 ) du 
and the source 5*o is given by 
q 2 (3-u 4 ) 



16u b 



1 



(A3) 



So 



° 2 4 

- g<7 u *05 



q 2 u 4 ft 



+ q 2 (l+U 4 )t qq 
32u 3 iqf 



32u 3 iq j 
3-u 4 



5o 



(A4) 



3" J 3-u 4 
Near the boundary Z Q has the following expansion 

■ 3 ' ^ - 4 ' ... (A5) 
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The zero mode of the boundary energy density is related 
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to Z®^ via [38 
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qR 



dq 



^ (f,0), (A7) 



with i?3 = | lim u _>o d^R and Rq = lim^^o R. 



3 The definition of Zq that we have given is valid only at oj = 
but this definition can be extended to nonzero ui. Upon doing so, 
the expression flAl) becomes much lengthier and, furthermore, 
one discovers that it features ratios of terms that both vanish at 
certain values of u between and 1. Although Zq is in principle 
well-behaved at these values of u, this feature makes Zq imprac- 
tical for numerical calculations when u/0. No such difficulties 
arise at uj = 0, which is the only case for which we use Zq. 
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